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Abstract The dynamical mass of a star cluster can 
be derived from the virial theorem, using the measured 
half-mass radius and line-of-sight velocity dispersion of 
the cluster. However, this dynamical mass may be a 
significant ovcrestimation of the cluster mass if the con- 
tribution of the binary orbital motion is not taken into 
account. In these proceedings we describe the mass 
ovcrestimation as a function of cluster properties and 
binary population properties, and briefly touch the is- 
sue of selection effects. We find that for clusters with 
a measured velocity dispersion of a\os ^ 10 kms~^ the 
presence of binaries does not affect the dynamical mass 
significantly. For clusters with ctios ^ 1 kms~^ (i-e., 
low-density clusters) , the contribution of binaries to ctios 
is significant, and may result in a major dynamical mass 
ovcrestimation. The presence of binaries may introduce 
a downward shift of Alog{Ly/Mdyn) — 0.05 — 0.4 in the 
log(Ly /Mdyn) vs. age diagram. 



using the equation derived by ISpitzeii (|l987l ): 
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1 Introduction 

An estimate for the mass of star clusters can be ob- 
tained from the virial theorem, using the projected half- 
mass radius i?hm and line-of-sight velocity dispersion 
(Tios- This dynamical mass estimate, Mdyn, is obtained 
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dyn ~ V 



G 



(1) 



where 77 « 9.75 is a dimensionless proportionality con- 
stant. The derivation of Afdyn using the expression 
above is valid under the following assumptions: (1) 
the cluster dynamics are described by the Plummer 
model, (2) all stars are single and of equal mass, (3) 
the cluster is in virial equilibrium, and (4) no selec- 
tion effects are present. Dynamical mass estimates 
for numerous cl usters have been obtained thi s way 
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Maraston et al.ll20 04: Bas tian et al.ll2006l : iLarsen et aL. 
20071 : IMoII et al.ll2007ll . When binary stars are present, 



however, Eq. ([T]) results in an ovcrestimation of the clus- 
ter mass. 

Observations have shown that the majority of the 
field stars are part of a b inary or multiple system 
( Duquennov fc Mavoil Il99ll IPischer fc Marcvl Il993 ). 



It is believed that most (if not all) binary stars are 



both observations ( 


Mathieu 1994: Mason et al. 


1998 


iKobulnickv & Frvei 


2007; iKouwenhoven et al.l 


2005 


i2007) and theory ( 


Goodwin & KrouDa 20051). 


Stars 



in binary systems exhibit not only motion in the gravi- 
tational potential of the cluster [particle motion CTpart), 
but also orbital motion Uoth in the binary system. 
Eq. ([1]) is applicable for ijpart (the centre-of-mass mo- 
tion of the binaries) , but results in an overestimation if 
(Tios, the superposition of cTpart and (Jorb, is measured. 
In this paper we therefore address the question: "How 
do binaries affect the dynamical mass of a star cluster?" 



2 Method and terminology 

We evaluate the effect of binarity on Mdyn using numer- 
ical simulations. We use the STARLAB package (e.g.. 
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Table 1 The default properties of the model used in our 
analysis, which we refer to as model R. In our analysis we 
vary the properties of modelled star clusters in order to 
find the effect of these changes on the derived dynamical 
mass Mdyn. At the bottom of the table we list the line-of- 
sight velocity dispersion aios of the individual stars, cTpart of 
the centre-of-mass motion of the binaries, and cTorb of solely 
the orbital motion of the binary components. Each value 
represents the width of the best-fitting Gaussian. 



Property 


Model R 


Model 


Plummcr 


Half-mass radius 


^hm = 5 pc 


Particles N = S + B 


iV = 18 600 


Total mass 


Mel = 10^ Mq 


Mass segregation 


No 


Virial equilibrium 


Yes 


Primary mass 


/Kroupa;0. 08-20 Mq{Mi) 


Binary fraction 


Fm = 100% 


Mass ratio 


f,{q) = 1; < g < 1 


Eccentricity 


fe{e) = 2e; < e < 1 


Orbital size 


/opik(a); 10 Rq- 0.02 PC 


Orientation 


Random 


aios (measured) 


1.20 kms^i 


Cpart (centres-of-mass) 


0.91 kms-i 


forb (binaries) 


0.14 kms^i 



When changing certain properties of the binary pop- 
ulation in a star cluster, such as the increasing the aver- 
age mass ratio (q) or binary fraction Fm , the total mass 
of the cluster increases. The latter increase results in 
larger particle velocity, which should not be fully at- 
tributed to the change in the distribution mass ratio 
distribution f(q) or in Fm. The increase in particle ve- 
locity is partially caused by the increased cluster mass. 
The latter affect can be compensated for by adjusting 
the number of particles N = S + B in the cluster. In 
this way the total cluster mass "before" and "after" the 
change in the binary population is equal. The effect of 
the changed binary property can now be studied accu- 
rately. Throughout these proceedings we keep the total 
mass for each comparison fixed to 10^ Mq, a typical 
value for open cluster-like objects. 

3 Dependence on cluster properties 

In the sections below we describe how the derived dy- 
namical mass depends on the structural parameters of 
a cluster and on its stellar mass distribution. 

3.1 Mass distribution and number of particles N 



Portegies Zwart et al.ll200ll ) to model star clusters with 
different structural and binary population properties. 
The default properties of our cluster model are listed 
in Tabledl we refer to these as 'model R'. We determine 
the best-fitting velocity dispersion from the velocity dis- 
tributions, and do not weigh by mass or luminosity. 

In the introduction we described three velocity dis- 
persions: the measured line-of-sight velocity dispersion 
(Tios, the particle (centre-of-mass) velocity dispersion 
Cpart, and the binary orbital velocity dispersion fiorb- 
Whether binarity is important depends on which of 
these (Tpait or tJorb dominate the measured fiios. We 
discriminate between three types of clusters: 

• Particle-dominated clusters. The measured ctios is 
dominated by cTpart; Eq. ([T]) is a reasonable approxi- 
mation and Afdyn/Mci < 110%. 

• Intermediate-type clusters. The contribution of dpart 
and CTorb to CTios is comparable; the dynamical mass 
overestimation is 110% < Mdyn/A^ci < 200%. 

• Binary-dominated clusters. The measured ctios is 
dominated by CTorb- The dynamical mass overesti- 
mation is significant: A/dyn/Afci ^ 200%. 

Typically, binary-dominated clusters have ctios ^ 1 
kms^"'^, while particle-dominated clusters have ctios ^ 
10 kms-i (see Sect.i]). 



Eq. 11]) assumes that all stars are single, equal-mass 
stars. In reality, the stellar mass spectrum is d efined 
by the mass function, such as the Kroupal ( 2001 ) IMF. 
The effect of the mass distribution on the dynamical 
mass determination, i n particular as a fun ction of time, 
is studied in detail bv iFleck et al.l (|2006r ). and will not 
be discussed further here. 

The number of particles N = S + B in a star cluster 
is related its total mass via M^i = {M)N, where (M) is 
the average mass of a particle (singles/binaries), which 
is defined by the mass distribution and the pairing prop- 
erties of the binaries. For a cluster consisting of single 
stars only, the dynamical mass given by Eq. ([T]); the 
determination of Afdyn is not dependent on the value 
of N. This is not the case if binaries are present. In 
this case, a larger N (or larger total mass) results in an 
increased CTpart , while the orbital motion of the binaries 
(reflected in CTorb) remains unaffected. Clusters with a 
large N are thus less affected by binarity, and allow for 
a more accurate determination of Afdyn than those with 
small N. 

3.2 The half-mass radius -Rhm and stellar density phm 

For a cluster with a given total mass and virial ra- 
tio, the velocity dispersion depends strongly on its size: 
fpart OC i?hm^^; scc Eq. ((T|) . Particlcs in a large (i.e.. 
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Fig. 1 The effect of the half-mass radius Rhm and density 
phm on the measured line-of-sight velocity dispersion ai^s 
and the overestimation of the dynamical mass of the cluster. 
The dashed lines indicate the results calculated from the 
centre-of-mass motion, and the solid curve the results for 
the measured line-of-sight velocities. 

sparse) cluster move slower in the potential than in a 
tight, high-density cluster. The orbital motion of bi- 
nary systems is independent of the cluster properties. 
In a low-density cluster the binary orbital motion is 
thus expected to dominate crios, while the presence of 
binaries is negligible in high-density clusters. The de- 
pendence on i?hm and the average density within the 
half-mass radius phm is illustrated in Fig. [TJ Note that 
the dynamical mass of a distant massive OB association 
(Mci « 10^ Mq] i?hm ~ 20 pc) can be overestimated 
by almost an order of magnitude, if its binaries are not 
properly taken into account. 

3.3 The virial ratio Q 




0.1 1 10 100 1000 0.1 1 10 100 1000 
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Fig. 2 The effect of the binary semi-major axis a on the 
measured line-of-sight velocity dispersion crios and the over- 
estimation of the dynamical mass of the cluster. All clus- 
ter properties are listed in Table [T] The solid curves show 
the results for clusters in which all binaries have an iden- 
tical a = ao, with ao along the horizontal axis. The dot- 
ted curve indicate the results for two realistic distributions: 
f(a) PC (Opik's law) and the IPuouennov fc Mavorl 

l|l99ll ) period distribution. The results for both distribu- 
tions are indistinguishable in this figure. The dashed lines 
indicate the results calculated from the centre-of-mass mo- 
tion. 

4 Dependence on binary population properties 

The reliability of a dynamical mass determination de- 
pends not only on the properties of the cluster (which 
are reflected in Cpart), but also the properties of the 
binary population (which are reflected in (Torb), in par- 
ticular on the size of the binary orbits and the binary 
fraction. The effect of each of these binary population 
properties is discussed in the sections below. 



After star formation the remaining gas from which the 
stars have formed is ejected by stellar winds and super- 
nova explosions. During this phase a star cluster may 
loose a significant fraction of its mass, resulting in a 
reduced gravitational potential. The cluster is now out 
of virial equilibrium, possibly even unbound, and starts 
expanding. Although Eq. ^ assumes that the cluster 
is in virial equilibrium, a correction for expanding or 
contracting clusters is easily made. Let Q = —Ek/Ep 
be the virial ratio, with Ek the total kinetic energy and 
Ep the total potential energy of the cluster. Clusters 
with Q < 0.5, Q = 0.5 and Q > 0.5 are contract- 
ing, in virial equilibrium, and expanding, respectively. 
For any value of Q, the relation between the true mass 
and the dynamical mass of a st ar cluste r is gi ven by 



Mci = (2Q)-iMdyn. iGoodwin fc Bastian. (|2006f ) define 



the effective star forming efficiency (eSFE) e as the star- 
forming efficiency that one would derive from the virial 
ratio under the assumption that the star-forming cloud 
was originally in virial equilibrium: Q = (2e)^^. Under 
this assumption, the dynamical mass overestimation is 
given by Mdy„/Mci = 



4.1 The semi-major axis distribution /(a) 

The distribution over semi-major axes (or periods) is 
one of the most important parameters that affect the 
interpretation of the observed uios, as the orbital ve- 
locity of a component is a binary system is propor- 
tional to a~^/^. Clusters containing wide binaries are 
less affected by binarity than those containing tight 
binaries. This effect is clearly shown in Fig. [H In 
binary-dominated clusters, which have ctios ~ CTorb, the 
velocity dispersion scales as ctios oc a^^/^, and there- 
fore Mdyn/Mc\ cx a^^. In the particle-dominated case 
(which occurs if most binaries are wide), the dynam- 
ical mass is a good representation of the true cluster 
mass. The intermediate case occurs approximately at 
log(a/i?hm) ~ ~4.7, where i?hm is the half-mass radius 
of the cluster. Note that the two most commonly used 
orbital size distributi ons, /(a) cx a^^ (Opik's law; e.g. 
Poveda fc AlleJ2 004 ) and the log-n ormal period distri- 



bution ( Duquennov fc MavoijllQQl ) practically give the 
same results. 
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Fig. 3 The effect of tfie binary fraction Fm on tiie mea- 
sured line-of-sigiit velocity dispersion crios and the overesti- 
mation of the dynamical mass of the cluster. Model param- 
eters are listed in Table [1] The dashed lines indicate the 
results calculated from the centre-of-mass motion, and the 
solid curve the results for the measured line-of-sight veloci- 
ties. 



4.2 The binary fraction Fj 



M 



The binary fraction is an important parameter, as it 
determines the relative weight that is given by Upart 
(the single stars) and dorb (the binaries) to the mea- 
sured value of (Tios • The dynamical mass of a star clus- 
ter with a low binary fraction is expected to be only 
barely overestimated. As the binary fraction increases, 
the overestimation gradually becomes larger. In the 
intermediate-case, the dynamical mass overestimation 
scales more or less linearly with Fm , which is illustrated 
in Fig. H 

4.3 The eccentricity distribution /(e) 

Stars in an eccentric orbit spend most of their time near 
apastron, where their velocity is small, and a small frac- 
tion of their time near periastron, where their velocity 
is large. The average velocity of a set of eccentric bi- 
naries at a certain point in time is therefore relatively 
small. For a cluster with highly eccentric orbits, the 
contribution of (Toib to uios is therefore smaller than for 
a similar cluster with circular binaries. The dynami- 
cal mass overestimation thus decreases with increasing 
average eccentricity. However, the effect of an incor- 
rectly adopted eccentricity distribution on the calcu- 
lated Mdyn is small, as compared to that of the un- 
certainty in /(a), Fm, and the selection effects. For 
model R our simulations show Mdyn/Afci — 80% for a 
cluster with circular binaries, while Mdyn/Afd = 50% 
if all binaries have e = 0.95. 



the massive primary star barely moves, while the low- 
mass companion orbits at high velocity. Depending on 
which star is measured, the orbital velocity contribution 
to a\os can be small or large. Most frequently the ve- 
locity of the brightest, and therefore most massive star 
is measured. A higher average mass ratio thus leads to 
a larger dynamical mass overestimation. However, in 
practice the effect of an uncertainty in the mass ratio 
distribution is much smaller than for an uncertainty in 
/(a), Fm and the selection effects. 



5 Selection effects 

Selection effects play an important role for the interpre- 
tation of the measured trios, in particular (i) the pro- 
jected radius from the cluster centre at which the mea- 
surement is performed and (ii) the stellar mass range 
that is included in the observations. The line-of-sight 
velocity dispersion of the centres-of-mass at a certain 
projected distance p from the cluster centre is 



par 



6414^ 



'1/2 



(2) 



( Heggie fc Hull 20031 ) . An expression for the dynamical 
mass overestimation as a function of p is obtained by 
substituting Eq. into Eq. H]): 



M, 



dyn 



M, 



V2 1 



i?2 



-1/2 



(3) 



For velocity dispersions measured in the cluster centre 
this results in a mass overestimation by ~ 40%. Mea- 
surements at the half-mass radius provide the correct 
Afdynj while measurements in the cluster outskirts re- 
sult in an underestimation of the mass. A possibly more 
important selection effect is introduced by the large 
brightness difference between stars of different masses. 
In reality, the determination of uios is dominated by 
the properties of stars in a certain mass range (more 
specifically, a certain brightness range). The measured 
fJios may therefore not be representative for the clus- 
ter as a whole. This may re sult in a further dynami- 
cal m ass overestimation (e.g., Kouwenhoven fc de Grijsl 
20071 1. A detailed analysis of the selection effects is 
necessary to properly take these selection effects into 
account. 



4.4 The mass ratio distribution f{q) 



6 When can binaries be ignored? 



Equal-mass stars in a binary system orbit each other 
at equal velocities. For a very low mass ratio binary, 



When studying a star cluster in detail, it is important to 
find out how reliable the measured dynamical mass is. 
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Fig. 4 The effect of binarity on the dynamical mass de- 
termination for a cluster with a mass Md and a half-mass 
radius Rhm- For young massive clusters, binarity generally 
has a small effect on Mdyn, while for distant OB associa- 
tions the overestimation of Mdyn can be up to an order of 
magnitude. 



Fig. m shows for clusters of different size i?hm and mass 
Mci whether they are particle-dominated, intermediate- 
case, or binary-dominated. The figure shows that 
the low-density OB associations are generally binary- 
dominated; their dynamical mass as obtained from inte- 
grated spectral lines would significantly overestimated 
their true mass. Young massive clusters, with masses of 
10^~^ M0 and typical half-mass radii of a few parsec, 
are of the intermediate case. Their dynamical masses 
are overestimated by a few per cent. From a practical 
point-of-view. Fig. [5] shows how the measured velocity 
dispersion uios should be interpreted, and how the re- 
sults depend on the intrinsic binary fraction. In the 
right-hand panel of Fig. [5] we additionally indicate the 
downwards shift Alog(Lv/Afdyn) = log ( Mdyn /Afci) in 
the luminosity-to-mass vs. age diagram as a results of 
binarity. 



7 Conclusions 

The total mass of a star cluster is often inferred from 
its line-of-sight velocity dispersion trios and half-mass 
radius i?hm, assuming virial equilibrium. The latter 
approach includes the assumption that no binaries are 
present. However, most stars are known to be in bi- 
nary systems, and their orbital motions provide an ad- 
ditional contribution to the measured crios- The latter 
value is now no longer representative for the (centre- 
of-mass) motion of the stars and binaries in the cluster 
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Fig. 5 Given the measured line-of-sight velocity dispersion 
(Jios, what is the effect of binarity? Can the effect of bina- 
rity on the dynamical mass derivation be ignored? And if 
not, how severe is the overestimation of the inferred dynam- 
ical mass? For a given trios, the curves indicate the differ- 
ence trioa — CTpart (left-hand panel), and the dynamical mass 
overestimation Mdyn/Md (right-hand panel). Results are 
shown for models with _Fm — 0% (dotted curves) , _Fm ~ 25% 
(dashed curves), and for 50%, 75% and fOO% (solid curves). 
We additionally indicate the shift A log(Lv/Mdyn) in the 
(Lv/Mdyn) vs. age diagram that is introduced by the pres- 
ence of binaries. 



potential. This effect may result in a significant dynam- 
ical mass overestimation. Depending on the magnitude 
of the dynamical mass overestimation, we distinguish 
between three types of clusters: particle-dominated 
(o-orb < CTpart), intcrmediate-case (trorb ~ Cpart), and 
binary-dominated clusters (aorb ^ fpart)- 

The orbital velocity of binaries is independent of the 
cluster properties (size, mass, etc.). Whether or not bi- 
nary motion affects trios is thus depends on the cluster 
properties. For clusters with a high stellar density (i.e, 
large Md or small i?hm), M^yn is generally unaffected. 
In the latter case trpait > 10 kms~^. The dynamical 
mass overestimation increases strongly with (i) a higher 
binary fraction and (ii) a smaller average orbital size. 
The dependence on the mass ratio distribution and ec- 
centricity distribution is small: A(Mdyn/Mci) < 5%. 
We additionally show that observing a certain subset 
of the cluster introduces a selection effects, which may 
result in a further mass overestimation by up to 40%. 
Analysis of the brightest (generally massive) stars in the 
cluster may further overestimate the dynamical mass. 
The dynamical mass overestimation can be negligible 
for the most massive clusters, while it may overestimate 
the true mass by up to an order of magnitude for sparse 
OB associations. The full results of this study have 
been presented in iKouwenhoven fc de Griid (|2007l ). A 
follow-up paper, which treats the selection effects prop- 
erly, is in preparation. 
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